4 th Fermi Symposium : Monterey, CA : 28 Oct- 2 Nov 2012 



1 



OVRO 40 m Blazar Monitoring Program: 

Location of the gamma-ray emission region in blazars by the study of 
correlated variability at radio and gamma-rays 

W. Max-Moerbeck 1 , J. L. Richards 2 , V. Pavlidou 3 , T. J. Pearson 1 , A. C. S. Readhead 1 on behalf of the 
Fermi LAT collaboration, T. Hovatta 1 , O. G. King 1 and R. Reeves 1 . 

1 California Institute of Technology, Pasadena, CA 91 125, USA 

2 Department of Physics, Purdue University, West Lafayette, IN 47907, USA and 

3 Max-Planck-lnstitut fur Radioastronomie, Aufdem Hugel 69, 53121 Bonn, Germany 



m 
o 



oo 

Oh 
i 

o 

in 



> 

m 
(N 

o 

m 

> 
X 



Blazars are powerful, variable emitters from radio to gamma-ray wavelengths. Even though the general picture 
of synchrotron emission at low energies and inverse Compton at the high energy end is well established, many 
important aspects of these remarkable objects are still not well understood. For example, even the location of 
the gamma-ray emission region is still not clearly established, with some theories locating it close to the black 
hole/accretion disk while others place it at parsec scales in the radio jet. Since mid-2007 we have carried out a 
large scale monitoring program at 15 GHz using the OVRO 40 m telescope. We are currently observing about 
1700 blazars twice per week. The sample includes all the Fermi-hAT detected blazars north of declination —20°. 
Here, we study the existence of correlated variability between these two bands for 86 sources bright enough to be 
detected weekly by LAT. The existence of correlated variability can be interpreted as an indication of a related 
spatial locations for the radio and gamma-ray emission, making the evaluation of its statistical significance a 
key goal of our program. A study of the statistical significance of these cross-correlations is presented along with 
a discussion of the Monte Carlo simulations used to evaluate them. More information about the conditions on 
the radio emission zone can be obtained through polarization monitoring which tells us about the configuration 
of the magnetic fields in this region. To study radio polarization variability we are building KuPol, a radio 
polarization receiver for the 12 to 18 GHz band that will replace the current total power receiver at the OVRO 
40 meter telescope. 



1. INTRODUCTION 

Blazars are a type of active galactic nuclei with a 
jet oriented close to the line of sight. Their spectral 
energy distribution is broad and characterized by the 
presence of two peaks: a low energy one produced by 
the synchrotron mechanism, and a high energy one 
produced by inverse Compton scattering of a photon 
field from the jet, accretion disk, the broad line region 
or the dust torus. Their emission is variable over a 
whole electromagnetic spectrum. 

Due to these characteristics, a complete under- 
standing of the emission mechanisms in blazars re- 
quires their study over a broad range of energies with 
simultaneous monitoring at different energy bands. 
This monitoring needs to be independent of the state 
of the source to avoid biases. Moreover, understand- 
ing the statistical properties of different classes re- 
quires the study of large samples of objects. 

Even though the general picture of synchrotron 
emission at low energies and inverse Compton at the 
high energy end is well established, many important 
aspects of these remarkable objects are still not well 
understood. Our work concerns one such question: 
where in the blazar structure is the gamma-ray emis- 
sion region located? Some theories locate it near 



the black hole/accretion disk Blandford & Levinson 



1995 , while others place parsecs away down the jet 



Jorstad et aL|2001 



directly image the gamma-ray emission regions in 
blazars, but we can investigate their location by study- 
ing the connection between the radio and gamma-ray 
variability. To enable such studies, we are carrying 
out a blazar monitoring program with a sample size 
and cadence matched to Fermi's capabilities. 

This paper presents a summary of the state of this 
project as was presented at the time of the sympo- 
sium. A detailed discussion concerning the data anal- 
ysis and results will be presented in Max-Moerbeck et 
al. (in preparation). 



2. OBSERVATIONS 

We are currently monitoring about 1700 blazars in- 
cluding all the sources north of —20° in declination 



in the Candidate Gamma-ray Blazar Sample Healey 



et al. 



with 



2008 



and all Fermi detected sources associated 
jlazars in the first and second LAT AGN catalogs 
Abdo et al. 2010a Ackermann et al. 20111 with the 



The Fermi Large Area Telescope (LAT) and other 
gamma-ray telescopes lack the resolution needed to 



same declination limit. The sources are observed in 
total flux density twice per week at a frequency of 15 
GHz with a 3 GHz bandwidth. A typical observation 
has about 4 mJy of thermal noise with an additional 
~ 3% uncertainty coming mostly from pointing errors. 
A complete description of the program can be found 
in Richards et al. 2011 . 

These observations complement the continuous 
monitoring of the gamma-ray sky in the band from 
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100 MeV to 300 GeV by Fermi-LAT. In this study, 
we consider all the gamma-ray sources detected in at 
least 75% of the monthly time bins as reported in 
Nolan et al. 2012 . This sample contains 86 sources 



for which we obtain weekly binned gamma-ray light 
curves for integrated photon flux in the energy range 
from 100 MeV to 200 GeV covering the time period 
from 4 August 2008 through 8 August 2011. The ra- 
dio light curves cover the period from 1 January 2008 
through 26 February 2012. 



3. MEASURING THE 
RADIO/GAMMA-RAY TIME LAGS 



To study the physical significance of cross- 
correlation peaks we use Monte Carlo simulations that 
account for the uneven sampling and noise properties 
of the light curves. The uneven sampling is incorpo- 
rated using the discrete correlation function Edelson 
fc Krolik|1988] and the noise properties of the sources 
are characterized by fitting their power spectral den- 
sity (PSD) with a simple power law in both bands 
(P(/) oc 1/ /^) using an adapt ation of the method 

2002] . Here we present re- 



presented in Uttley et al 



suits of a study using populations values of the PSD 
power law index at both bands. For the gamma-ray 
band we adopt the average of the results obtained for 
the mean PSD of the brightest BL Lacs and FSRQs 
(/3 7 = 1.6 from |Abdo et al-1|2010b]). For the radio 



band we use the mean population value we obtain in 
this work from fits to the radio light curves of the 
complete OVRO monitored sample (/3 ra dio = 2.3). 

The significance of a cross-correlation peak is mea- 
sured by obtaining an estimate of the distribution 
of random cross-correlation amplitudes at the given 
time lag. This distribution is estimated by simulat- 
ing a large number of independent time series (we use 
20,000 in this case) with the appropriate power spec- 
tral densities and observational noise. 



5. THE CASE OF MRK 421 

The well known blazar Mrk 4 21 underwent a major 
gamma-ray flare in July 2012 D'Ammando & Ori- 



|enti|2012| which was followed by a major event at 15 
GHz, the largest flare since at least 1980 as shown by 
a comparison to University of Michigan Radio Astron- 
omy Observatory 14.5 GHz data Hovatta et al.pOF? 



When we extend the light curves to include the flar- 
ing period we are able to constrain the power spectral 
density in both bands and find correlated variability 
with > 2cr significance. Preliminary light curves and 
results of the significance analysis including only the 
rising part of the radio flare are shown in Figure [3] 
This result has to be considered with care as in this 
case we are extending the period analyzed after ob- 
serving the flare, thus doing "a posteriori" statistics. 



6. KUPOL: RADIO POLARIZATION 
MONITORING 

Blazars produce their radio emission by the syn- 
chrotron mechanism which can result in polarized 
emission. Variations in the polarization can help us 
understand the evolution of the magnetic fields in the 
emission regions and its relation to high-energy emis- 
sion. To study these magnetic fields we are building a 
new receiver which will also measure linear polariza- 
tion in the Ku-band. This new receiver, which we call 
KuPol, will cover the band from 12 GHz to 18 GHz 
with 16 MHz spectral resolution, which can be used 
for examining spectral indices and removal of radio 
frequency interference. This will be an important ad- 
dition to our program and commissioning is planned 
for spring 2013. 



7. SUMMARY 



4. SIGNIFICANT CORRELATIONS 

Of the 86 sources with good cadence in radio and 
gamma-ray bands included in this study, only 63 show 
variability in both bands and are thus included in the 
cross-correlation analysis. Among these 63 we find one 
case with > 3a significant correlations, J0238+1636 
(Figure [T] 0.2 cases are expected by chance). There 
are 6 other cases with > 2a significant correlations 
(3 cases are expected by chance). However, we con- 
sider five of these unreliable because of high noise or 
slow trends. The remaining example, J1504+1029, is 
shown in Figure [2j 



We study the connection between radio and 
gamma-ray emission in Fermi detected blazars and 
find that significant correlations are present in a 
small number of these sources with the current data 
set. The major flare on Mrk 421 shows correlated 
radio/gamma-ray activity. These cases add support 
to the idea that radio and gamma-rays are produced 
in related regions for some sources, but more data are 
needed to test whether this holds for the entire pop- 
ulation. In particular longer radio and gamma-ray 
light curves will allow us to probe the variability on 
different time scales, and to incorporate information 
for multiple flaring periods. This will improve the 
statistics and the significance of possible correlations. 
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Figure 1: Light curves and cross-correlation results for J0238+1636. The most significant peak it at —30 ± 8 days and 
has a 99.9% significance. Color contours on right panel indicate the significance of the cross-correlations with red for 
la, orange for 2a and green for 3a. 
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Figure 2: Light curves and cross-correlation results for J1504+1029. The most significant peak it at — 40 ± 13 days and 
has a 98.5% significance. Color contours as in Figure [I] 
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Figure 3: Light curves and cross-correlation results for Mrk 421. There are two significant peaks one at 520 days with 
98.8% significance and a shorter one at —40 days with 96.3% significance. Color contours as in Figure [l] This lag 
determination needs to be confirmed with more data including the decaying part of the radio flare. 
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